International 

open 9 access Journal 

Of Modern Engineering Research (IJMER) 

I ~1 



Magneto hydrodynamic Rayleigh Problem with Hall Effect and 
Rotation in the Presence of Heat Transfer 

P T Hemamalini 1 and M Deivanayaki 2 

1 Professor , Department of Science and Humanities , Faculty of Engineering, Karpagam University, Coimbatore, 
2 Research Scholar, Department of Mathematics, Karpagam University, Coimbatore. 



\ 

ABSTRACT: The main objective of the paper is to investigate the magneto hydrodynamic flow version of 

the Rayleigh problem including Hall effect and Rotation in the presence of Heat transfer. Analytical 
solution has been found depending on the physical parameters including the Hall parameter N , Hartmann 

number M f Grashof number Gr ? Prandtl number Pr and the Rotation parameter K 2 . The influence of 
these parameters on velocity profiles and temperature profiles are demonstrated graphically and the results 
are discussed. Further, it is observed that increase in the Hall parameter and Rotation parameter leads to 
decrease in the velocity profiles. It is also found that the increase of Grashof number results in the decrease 
of primary flow and increase of secondary flows. 
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I. INTRODUCTION 

In recent years, the analysis of hydromagnetic flow has applications in diverse fields of Science and 
Technology such as soil sciences, astrophysics, nuclear power reactors etc. The study of MHD flow problems 
has achieved remarkable interest due to its application in MHD generators, MHD pumps and MHD flow meters 
etc. Geophysics encounters MHD phenomena in interaction on conducting fluids and magnetic fluids. The 
rotating flow of an electrically conducting fluid in presence of magnetic field has got its importance in 
Geophysical problems. The study of rotating flow problems is also important in the solar physics dealing with 
the sunspot development, the solar cycle and the structure of rotating magnetic stars. 

Stoke’ s analysed the flow of an incompressible viscous fluid past an impulsively started infinite 
horizontal plate in its own plane. Rossow[l] initiated Rayleigh’s problem for non-conducting plate while 
Chang&Yen[2] have taken the plate to be perfectly conducting in a various transverse magnetic field. The effect 
of coriolis forces on Rayleigh Problem was considered by Sathi [3]. Rayleigh problem in MHD with suction 
was considered by Girish Chandra Pandey [4] . 

The Hall effect is due merely to the sideways magnetic force on the drifting free charges. The electric 
field has to have a component transverse to the direction of the current density to balance this force. In many 
works on plasma physics, the Hall effect cannot be ignored as it has a significant effect on the flow pattern of an 
ionized gas. Hall effect results in a development of an additional potential difference between opposite surfaces 
of a conductor for which a current is induced perpendicular to both the electric and magnetic field. This current 
is termed as Hall current. The effect of Hall current on MHD Rayleigh’s problem in ionized gas was discussed 
by Mohanty [5]. Solution of Rayleigh problem for conducting fluid was considered by Abd-el-Malek etal[6]. 
Magneto hydrodynamic Rayleigh problem with Hall effect was studied by Haytham Sulieman[7]. Deivanayaki 
etal[8] studied the MHD Rayleigh problem with Hall effect and Rotation. In this study the effect of the Hall 
current and rotation on the magneto hydrodynamic flow version of the classical Rayleigh problem in the 
presence of Heat transfer was considered. 

II. FORMULATION OF THE PROBLEM 

The scenario under investigation comprises an incompressible electrically conducting, viscous fluid 
past an infinite vertical plate occupying the plane y = 0. The x-axis is taken in the direction of the motion of the 
plate and z - axis lying on the plate normal to both x and y - axis. Initially it is assumed that the plate and the 
fluid rotate in unison with a uniform angular velocity Q about the y - axis normal to the plane are at the same 
temperature T everywhere in the fluid. At time t > 0, the plate starts moving impulsively with the uniform 
velocity in its own plane along the x-axis. Also the temperature of the plate is raised/lowered to T^. A uniform 
magnetic field H 0 , parallel to y - axis is imposed. The schematic diagram is given by 
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Fig 1: Schematic representation 



Under the above assumptions, in the absence of an input electric field the governing boundary layer equations 



are 

Equations of continuity V. q = 0 (2.1) 

Equation of motion 

^ + (q. V)q + 2£2 x q = — - VP + yV 2 q + - J x H + pg (2.2) 

The energy equation 

KV 2 T = pC v ^- t (2.3) 

The generalized Ohm’s law, neglecting ion-slip effect but taking Hall current is, 

J -= (E + qxH)~— (2.4) 

c n.e 



6^ TH 

where a = (is the electrical conductivity). Here J is the current density, t is the time, p is density, y is 



kinematic viscosity, e is electric charge, m is mass of an electron, n is the electron number density, t is the mean 
collision time and p is magnetic permeability. 

Let T - T x = 9(y, t) 

The initial and boundary conditions are 
t < 0: u = 0, v = 0, 0 = 0 for y > 0 
t > 0:u = U 0 ,v = 0, 6(0, t) — ae l0)t fory = 0 
u = 0,v = 0, 0(oo, t) = 0 as y -> oo 

H x -> 0 ,H y = H 0 , H z -> 0 as y -> oo (2.5) 

At infinity the magnetic induction is uniform with components (0,//,0), and hence the current density vanishes 
and since the free stream is at rest, it follows from generalized Ohm’s law that E = 0 as y -> oo. Assuming small 
magnetic Reynolds number for the flow, the induced magnetic field is neglected in comparison to the applied 
constant field H 0 . 

Introducing the non-dimensional quantities: 



u o-y „ ,* “ „* * wqi Q . k ^ r 

y = , u = — , v = — , t = — ,6 = -,Gr = 



ukt 



U 0 



U o' 



gpyq 



,N = cot ,M 2 = 



pus 



Pr = r^ K 2 = r^ r» 

*■ put Y u o 2 

All the physical variables are defined in the Nomenclature. 

Equations (2.1), (2.2) and (2.3) transform to the following non-dimensional forms, respectively. 



du 



d z u 



_ _ _ _ r T n _ 1V1 

dt dy 2 ^ (i+^2) 

dv d 2 v M 2 

dt ~ dy 2 + (1+V 2 ) 
d 2 6 _ d6_ 

dy 2 ^ dt 



(u + Nv ) — 2vK 2 



(Nu - v) + 2uK 2 



The corresponding boundary conditions in non-dimensional form are 
For t < 0: u(y, 0) = v(y, 0) = 0, 6(y, 0) = 0 for all y 
For t > 0: u(0, t) = 1, v(0, t) = 0, 6(0, t) = e iy/t for y = 0. 
u(y, t) = 0, v(y, t) = 0, 6(oo, t) = 0 as y -> oo 



( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 



( 2 . 10 ) 
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III. METHOD OF SOLUTIONS 



The equations (2.7) and (2.8) can be combined using the complex variable q = u + iv 

57 = 0-[(^( 1 - iw >- 2iK V Cri ' 

The boundary condition (2.10) are transformed to 
q(y,0)=0, q(0,t) = 1, q(co,t) = 0 

Q(y, 0) = 0, 0(O,t)= 0(°o,t) = O 

Substitute 6{y, t) = e n//t f (y) in equation (2.9), we have 

/"GO - iysPrfiy) = 0 

The equation (3.3) can be solved under the boundary condition 

/( 0) = 1, / (°o) = o 

Hence the solution is f(y) = e^ y ^ u/,Pr 

Separating the real and imaginary parts and taking the real part only we get 
6(y,t) = e~ yxl Cos{yA - x x y) 




Also substitute q(y, t) = e iy/t g(y ) in (3.1), we have 

g"iy) - (.w+P)g(y) = -Gre-yJ 1 ^ 

The equation (3.7) can be solved under the boundary conditions, 
0 (O) = e-‘*, 0 (oo) = O 
Therefore the solution is 



(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 



g(y) = ( 



+ . 



Gr 



o-y^Jiys+P 



Gr 



i-yfiy/Pr 



iy/Pr-(yA.+p)J iy/Pr-(iy/+p ) 

Separating equation (3.9) into real and imaginary parts, we have 



u — e 



-yx 3 



Cos(y4:)Cos(ifd — yx 4 ) + e yX3 Sin(igt)Sin(y4: — yx 4 ) + 



(y/Pr-y/-b)Gr e yx ^Sin {yA.-yx\) 



a Gr e yx ^Cos{yA:-yx/i t ) (y/Pr-y/-b)Gr e yxi Sin(\yt-yxi) a Gr e yxi Cos(iyt-yx i) 



( y/Pr — y/—b) 2 +a 2 



( y/Pr — y/—b) 2 +a 2 



( y/Pr — y/—b) 2 +a 2 



+ ■ 



( y/Pr — y/—b) 2 +a 2 



v = e yx 3 Cos(igt)Sin(igt — yx 4 ) — e yX3 Sin(igt)Cos(igt — yx 4 ) — 



(y/Pr-y/-b)Gr e yX 3 Cos(yrt-yx4) 



a Gr e yx ^Sin {yrt— yx 4 ) , ( y/Pr — y/—b)Gre yxi Cos(yrt—yx i) , a Gr e yxi Sin(yrt—yx i) 

I , „ , N T . 1 I 



( y/Pr — y/—b) 2 +a 2 



(y/Pr — y/— b) 2 +a 2 



( y/Pr — y/—b) 2 +a 2 



(y/Pr — y/— b) 2 +a 2 



(3.9) 



(3.10) 

(3.11) 



The constants involved in the above discussion have been obtained but not presented here for the sake of 
brevity. 

In order to attain a physical insight into the problem, we have carried out numerical calculations for the velocity 
fields and temperature field at the plate due to time, Hartmann number, Hall parameter, Rotation parameter, 
Prandtl 
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FigS: Variation of u when Pr increases with y Fig9: Variation of v when Pr increases with v 
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FiglO: Variation of u when rotation increases with y Figll: Variation of v when rotation increases with y 




IV. CONCLUSION 

Fig. 2, 4 and 6 depict the variation of the velocity component u against y under the influence of 
Hartmann number, Hall parameter and Grashof number of heat transfer respectively. The velocity component 
decreases in all the cases. From fig. 3 and 7 it is clear that the velocity component v increases against y under 
the influence of Hartmann number and Grashof number. In fig. 5, the effect of Hall parameter on the velocity 
component v against y are shown and v decreases. 

We see from fig. 8 and 9 that u increase and v decrease when the effect of Prandtl number is 
considered. Fig. 10 and 11 shows that the velocity components u and v decreases with the increase of Rotation 
parameter against y. The temperature profile decreases with the increase of time against y as shown in Fig. 12. 
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